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SUMMARY 

The  problem,  investigated  is  the  resclvability,  in  the  presence  of  noise, 
of  point  sources  vhich  are  separated  in  azimuth  by  less  than  the  width  of 
the  main  lobe  of  the  gain  pattern. 

The  probability  of  correctly  resolvin^  two  sources  is  derived  for-  a 
particular  class  of  decision  methods  as  a  function  of  the  strengths  of  the 
sources,  their  angular  separation  and  noise  level.  Numerical  results  are 
presented  for  a  case  vhich  corresponds  more  accurately  to  optical  or  infra¬ 
red  uevices  than  to  radar. 
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I .  INTRODUCTION 


The  question  of  resolvability  of  point  targets,  or  point  sources,  arises 
in  radar,  infrared,  and  optical  applications. 

First  let  us  distinguish  the  problem  of  resolution  from  that  of  accuracy. 

The  accuracy  problem  is  concerned  with  the  errors  in  determining  the  position 

(l  2 ) 

of  an  isolated  source.  It  is,  for  example,  well  known  ’  that  the  error 
in  determining  angular  position  of  a  radar  target  may  be  much  less  than  the 
width  of  the  main  lobe  of  the  antenna  beam. 

The  word  "resolution"  covers  several  related  concepts.  In  photography 
or  radar  ground  mapping,  resolution  has  to  do  with  the  quality  in  some  sense 
of  the  picture,  or  with  the  ability  to  pick  out  certain  features  of  the 
picture.  On  the  other  hand,  the  type  of  resolution  with  which  we  shall 
here  be  concerned  is  the  following:  suppose  the  received  signal  is  known 
to  be  due  to  a  relatively  small  number  of  point  sources;  resolution  consists 
in  deciding  exactly  how  many  sources  are  actually  there.  We  will  also  limit 
the  discussion  to  angular  resolution;  resolution  in  range  or  range  rate  will 
not  be  considered, 

A  widely  used  rule  of  thumb  is  that  point  sources  are  resolvable  if 
they  are  separated  by  (roughly)  at  least  the  width  of  the  main  lobe  of  the 
gain  pattern.  (By  gain  pattern  is  meant,  in  radar  applications,  the  antenna 
gain  pattern;  in  infrared  or  optical  applications,  the  intensity  pattern  of 
the  image  of  a  point  source.  )  While  this  is  a  useful  rule  of  thumb,  it 
falls  to  answer  many  important  questions.  This  can  be  illustrated  by  the 
fact  that,  in  the  absence  of  all  noise,  and  if  the  gain  pattern  were  per¬ 
fectly  known,  it  would  in  many  cases  be  possible  to  make  error-free  decisions 
as  to  the  number  of  sources  regardless  of  their  relative  positions  or  strengths. 
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A  re  solvability  criterion  should  therefore  involve  the  probabilities 
of  making  correct  decisions  as  to  the  number  of  sources  present,  and  the 
manner  in  which  these  probabilities  depend  on  the  noise  level,  on  the  strengths 

and  positions  of  the  sources,  and  on  the  gain  pattern. 

6 

Ideally,  the  criterion  should  reflect  the  intrinsic  dependence  of  re- 
solvability  on  these  factors,  where  by  "intrinsic"  is  meant  the  (in  some 
sense)  irreducible  dependence,  apart  from  the  particular  decision  method 
employed.  However,  due  to  difficulties  in  arriving  at  a  satisfactory  treat¬ 
ment  of  "intrinsic"  resolvability,  the  approach  to  be  used  here  is  based  on 
analysis  of  the  performance  of  a  particular  class  of  decision  methods. 

It  is  clear  that  the  resolution  problem  as  defined  above  i3  a  problem 
of  statistical  hypothesis  testing.  However,  the  problem  of  finding  an 
optimum  test  cannot  be  answered  by  she  Neynan- Pear son  likelihood  ratio  test, 
since,  if  the  problem  is  not  to  be  simplified  out  of  existence,  the  positions 
and  strengths  of  the  sources  must  be  regarded  aa  not  known  a  priori.  Thus, 
the  problem  is  one  of  testing  composite  hypotheses  against  composite  alter¬ 
natives,  leading  to  all  the  well-known  difficulties  of  defining  optimum 
tests  in  “situations  of  this  type. 


II.  A  CLASS  OF  RESOLUTION  DECISION  METHODS 

It  will  be  assumed  that  we  are  dealing  with  a  detection  device  scanning 
in  azimuth.  The  received  signal  is  assumed  to  be  due  to  at  most  two  point 
sources.  It  will  be  assumed  that  the  output  of  the  detection  device  is 
sampled  at  a  aeries  of  angular  positions  x^,  i  =»  1,  •••,  K,  resulting  in  the 
presentation  to  a  decision-making  device  of  a  series  of  magnitudes  v.,  i  =>  1, 
•••,  K.  (If  we  are  dealing  wit):  a  pulsed  radar,  this  sampling  occurs  auto¬ 
matically;  otherwise,  such  a  sampling  may  always  be  defined  to  be  part  of 
the  decision  procedures  which  are  to  be  considered.  ) 

It  will  be  supposed  that,  when  two  sources  are  present,  one  can  repre¬ 
sent  by 

v.  >  ^  F(x.  -  tl)  ♦  a,  F(Xl  -  t2)  h  nr  (l  ) 

Here  and  are  equivalent  to  the  strength  of  sources  located  at  angular 

positions  and  t^;  F(x)  represents  the  gain  pattern;  and  n^  is  the  noise 

component  of  v  .  The  single-source  case  may  be  obtained  by  putting  either 

0^  *  0,  ^  °  0  orr^aXg-  (Not  all  physical  situations  can  be  represented 

by  Eq.  1;  more  will  be  said  about  this  below.  ) 

0 

The  following  assumptions  will  be  unde  regarding  F  and  n.  : 

(a)  F  is  a  known  function.  Thus,  attention  is  focussed  on  tre  limitations 
on  resolvability  due  to  noise;  that  due  to  imperfect  knowledge  of  the  gain 
pattern  is  ignored. 

(b )  The  gain  pattern  F  is  assumed  to  be  a  non-negative,  even  function. 

(c)  The  range  of  observations  v  ,  i  =  1,  •••,  K,  the  gain  pattern  F,  the 
possible  source  positions  t  and  the  density  of  observations  are  such  that, 


to  good  approximation 
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(d)  The  noise  components  n^  have  known  means  n.  (hpnce,  the  means  are  inde¬ 
pendent  of  the  unknown  parameters  a^,  ci0,  t^,  t^).  Thus,  it  is  no  loss  of 
generality  to  assume  that 


n.  B  0,  all  i 

i  s 
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It  will  not,  however,  be  assumed  that  the  higher  order  moments  of  n,  are 
necessarily  independent  of  a  ,  a.0,  t^,  t0. 

(e)  The  noises  n.  and  n,  are  statistically  independent  for  i  ^  J. 

These  assumptions,  incidentally,  simplify  many  of  the  mathematical  ex¬ 
pressions  but  are  not  really  essential  to  the  basic  approach. 

Now,  define 
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Observe  that,  if  =  0,  all  i,  then 
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A  reasonable  decision  method  vould  therefore  be  to  announce  two  sources 
or  one  source  according  to  whether  V  is  greater  or  less  than  a  pre -assigned 
threshold.  Observe  that,  if  n  ®  0,  all  i,  then  V  is  zero  if,  and  only  if, 
either  •»  0,  a  0,  or  ■>  t  . 

It  turns  out,  however,  that  this  decision  method  must  be  modified  in 
most  cases  of  interest,  for  the  following  reason:  if  n  ^  0,  then  (denoting 
expected  values  by  bars  ) 
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In  some  interesting  cases,  the  quantities  rm  depend  on  and  a0. 

Also,  the  higher  moments  of  V  as  defined  by  Eq.  (7)  may  depend  on  ci^  and  ou, 
even  if  the  statistics  of  n^  do  not. 

The  class  of  decision  methods  to  be  analyzed  can  be  defined  as  follows: 
announcement  of  two  sources  or  one  source  according  to  whether  V  is  greater 
to  or  less  than  a  pre-assigned  threshold,  where  V  is  some  function  of  and 
Vr  .  The  specific  form  of  V  will  depend  on  the  manner  in  which  the  moments 
of  n,  depend  on  ol,  a,  t, ,  and  t„,  in  the  particular  case  under  consideration. 

1  1  2  1  c 

It  is  clear  that  all  necessary  information  about  the  performance  of 
this  class  of  decision  procedures  can  be  deduced  from  the  joint  probability 
density  of  and  V^. 
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III.  THE  JOINT  PROBABILITY  DENSITY  OF  V  AND  V, 


It  will  now  be  assumed  that  the  number  K  of  observations  is  sufficiently 

great  so  that  both  and  are  normal  random  variables  (within  the  range 

# 

of  probabilities  which  we  wish  to  compute).  Thus,  the  joint  density  of 

-  -  2  2 

and  VQ  will  be  completely  specified  if  one  cam  compute  V^,  V^,  , 

and  V^V  ,  or  the  equivalent.  Computation  of  vhese  quantities  is  a  straight¬ 
forward  process.  The  results  are, 
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and,  writing  for  convenience 
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This,  of  course,  does  not  imply  that  V  is  normal. 
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IV.  RESULTS  FOR  A  PARTICULAR  CASE 


In  order  to  proceed  further  it  is  necessary  to  make  definite  assumptions 
a^  to  the  moments  of  n  ,  which  depend  on  the  particular  detection  device 
under  consideration. 

As  an  illustrative  case,  it  will  be  assumed  that 
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This  (with  suitable  normalization  )  would  be  the  case  if  the  n^  were  Gaussian 
with  means  and  standard  deviations  independent  of  a^,  a0.  t^,  and  x^.  This 
might  describe  some  infrared  or  optical  detection  devices;  it  is  not  so  easy 
to  + link  of  radar  cases  which  would  be  adequately  described  in  this  manner. 

The  first  and  second  order  moments  of  and  for  this  case  can  be 
obtained  by  inserting  (lp)  into  (9)  -  (l^-);  the  results  need  not  be  written 
out  in  detail.  It  turns  out  that 
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An  appropriate  resolution  test  is  defined  by  V  >  threshold,  where 
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Note  that  in  the  absence  of  noise,  one  would  have 
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Even  with  all  the  special  assumptions  which  have  thus  far  been  made, 
it  would  still  be  possible  to  generate  numerical  results  ad  infinitum. 
Therefore,  numerical  results  will  be  shown  only  for  the  following  special 
case : 


Let 


2 

F(x)  -  e"X 


(19) 


(20) 


Then,  approximately, 


B  -  20  v/T" 


(21) 


(22) 


The  range  of  observation  will  be  chosen  from  x  -  -2  to  x  -  +2.  With 
this  choice  of  observation  range,  K  ^  50.  The  choice  of  observation  range 
ought  really  to  be  the  solution  of  an  optimization  problem.  However,  the 
choice  is  here  somewhat  restricted  by  the  requirement  that  the  sums  in  Eq's. 
(2)  and  (3)  should  satisfy  Eq's.  (2)  and  (o)  to  good  approximation,  for  all 
values  of  and  under  consideration.  We  will  be  interested  in  values 
of  -  t2  up  to  about  1.5,  and  the  range  x  -  -2  to  x  ■  +2  was  chosen  as, 
roughly,  the  smallest  range  for  which  (2)  and  (3)  would  hold  for  all  values 

of  -  t2  up  to  1.5. 

The  results  obtained  foi  this  case  will  cleanly  also  apply  to  any  case 


of  the  form: 


F(x  )  =  ex:  -  my  I  ;  observation  range  :<  «  ~2\‘  to  .<  *  +  <rf';  x  -  » 

\  1+1 

and  all  i  values  rr.ultipl.ed  by  y. 

/II  parut.eter  value;,  will  re  chosen  sc  shut  t.  e  probability  of  failure 
to  let.eot  any  source  nt  all  la  negl bio . 

Figures  L-10  show  prot  ability  d.str:  :  .it  Ion  fur.c  1 1  o:  .3  for  V  (as  del-:.'. el 
by  (l/))  for  a  variety  of  val  :es  of  u  ,  a.  ,  and  1  -  ~  .  In  these  fi, yores. 
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sources,  if  in  reaxity  only  one  source  is 
present  character1’ zed  by  *he  parameter 
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Frci;.  Fig.  1,  we  •  ee  tint  •  r.  the  .resent  co.ro .  p[  two  1  '1,  ,  C]  is  nearly  independent 
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of  1  ,  provided  t:.  s  j.  r  rabil.tr  :s  relatively  sr.-all. 

Figures  11,  12,  o:.d  1.'  sr.c-w  P[  two  i  a,  u  ,  i,  -  1  ]  subject,  to  the  con- 
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It  ran  also  be  verified  :  the  particular  cane  under  consideration  that 
elution  j  er fonrar.ce ,  when  ;  resented  in  the  for::,  of  curves  like  those  of 
s.  *.1-1'',  re:. '.a .ns  roughly  tie  sat.e  if  x  ,  -  x.,  Q  and  rx-!~  a"e  all 
tijlied  by  ti.e  sn...e  :  actor  (otr.er  u^rarvters  reraaininp  unaltered). 
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v_1__  ADDITIONAL  TOPICS  PGR  I  INVESTIGATION 

Coiae  additional  volke  which  would  *..».•  usel.il  tc  investigate  are  as 
follows : 

i.  Investigation  cf  cure  ;  l.i  w..,  4  :.e  ’.or  ,>ls  cf  n.  (ot  ,er  than  the 

first)  depend  on  Clj  ,  O^,  t  ,  and  t„.  Th  is  would  include  some  radar  cases; 
for  example: 

(a)  Noise -like  targets  (or  noise  'ai.jr.ing )  detected  by  a  radar. 

In  ti.is  case,  the  output  cf  a  square  law  envelope  detector,  properly  norm!- 
ized,  would  have  probability  density  function 
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(b  )  Non -fluctuating;  targets,  tire  re  la'  ive  phase  of  one  target 
return  with  respect  to  the  ether  varying  .'audomiy  fret;,  pulse  to  pulse. 

One  would  expect  that  in  these  cases  the  resolvability  of  targets 
separated  by  less  than  a  beuinwidth  would  be  seriously  degraded  as  compared 
with  the  case  considered  in  the  previous  section. 

2.  Dropping  the  requirement  that  F  bq  an  even  function  and  always 
positive.  This  might  be  of  interest  in  analyzing  a  scanning  monopulse  radar 
Here  it  might  be  necessary  to  use  functions  of  the  data  ether  than  arid 
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